Comment on "Kondo resonances and Fano antiresonances in transport through 

quantum dots" 
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We show that the numerical method [M. E. Torio, K. Hallberg, A. H. Ceccatto, and C. R. Proetto, 
Phys. Rev. B 65, 085302 (2002)] does not reproduce correctly the Kondo physics in quantum dot 
systems with interaction. 



In a recent paper, M. E. Torio, K. Hallberg, A. H. Cec- 
catto, and C. R. Proetto (THCP) presented a numerical 
calculation of the conductance of an interacting system 
consisting of a quantum dot coupled to the leadsi Their 
results exhibit a fair qualitative agreement with the cor- 
responding experiments by J. Gores et al& However, a 
closer quantitative inspection of the results concerning 
the Kondo regime clearly shows that their method fails 
to reproduce adequately the main features of the conduc- 
tance in this regime. 

Numerical method of THCP consists of exact diago- 
nalization of a small cluster containing the region with 
interaction, coupling the cluster to external reservoirs 
(leads) , and finally self-consistent calculation of the corre- 
sponding Dyson equation. As pointed out in THCP, the 
method should be reliable as long as the corresponding 
Kondo cloud does not exceed the size of the cluster. In or- 
der to meet these criteria, THCP chose t' ' jt = 1/V2 and 
estimated (without giving arguments) the correspond- 
ing Kondo correlation length of about 10 lattice sites for 
U/A = 12 with A = 2t l2 /t. Using this set of parame- 
ters for the Anderson model^ they calculated the conduc- 
tance of substitutional quantum dot configuration. For 
the symmetric situation, eo = —U/2, their conductance 
reaches the unitary limit, G = Gq = 2e 2 /h. They fur- 
ther claimed that this is a nontrivial check of the method 
which proves that their finite-system approach produces 
a Kondo peak with the exact spectral weight at the Fermi 
energy^ 

In this Comment we show that although the conduc- 
tance of THCP reaches the unitary limit, the method 
does not reproduce correctly the Kondo peak. In 
Fig. 1(a) the conductance of THCP is shown compared to 
the result obtained with the two-point sine conductance 
formula from Ref.3using the variational wave function of 
Gunnarsson and Schonhammerai together with the con- 
strained path quantum Monte Carlo method (CPMC)i£*£ 
In the empty orbital regime, |eo + U/2\ > U/2, all results 
match perfectly. However, for |eo + U/2\ < U/2, except 
in the symmetric point, the result of THCP drastically 
fails to reproduce the main fingerprint of Kondo physics, 
i.e., the plateau in the conductance in the Kondo regime. 
The reason is the low energy resolution which in turn is 
a fingerprint of finite-size effects of a small clusters used 
by the method of THCP. It should be pointed out that 



even relatively large i! cannot remedy this deficiency of 
the method, because the problem falls into the unusu- 
ally narrow bandwidth regime 2t <C U, and e Q <C — 2t, 
eo + U >• 2t, where the scattering processes which in- 
volve electrons or holes in the band edges cannot cause 
real charge fluctuations of the impurity^ resulting in the 
reduction of the Kondo temperature. On the one hand, 
narrow bandwidth should improve the energy resolution 
of the method^ on the other hand, however, the Kondo 
temperature T% in this limit is reduced^ 
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Smaller t (or larger t') thus cannot bring the method 
of THCP into the regime of reliability unless the num- 
ber of sites in the cluster is significantly increased or the 
strength of the interaction, U/A, is reduced. 

In a more common limit of a wide bandwidth, 2t 3> 
L7(i& the conductance can be calculated exactly using 
the Friedel sum rule and the Bethe ansatz result pii In 
Fig. 1(b) the exact conductance for wide bandwidth is 
presented, together with the result of the two-point sine 
formula, using variational wavefunction and the CPMC 
method. The results agree for all regimes. Small discrep- 
ancies can be attributed to the finite bandwidth of both 
numerical methods (t = 1.04E7) and to the variational na- 
ture of the wave functions used. We would like to stress 
that due to the lower Kondo temperature in the narrow 
bandwidth regime, Fig. 1(a), the Kondo plateau is obvi- 
ously more pronounced compared to the wide bandwidth 
limit, Fig. 1(b). 

In Fig. 1(b) the conductance using an elementary 
Hartree-Fock calculation is also presented. It is well 
known that this approximation cannot adequately repro- 
duce the Kondo phenomenon. Nevertheless, the unitary 
limit in the symmetric point, eo + U/2 = 0, is in this 
approach as well reproduced. It is a consequence of the 
particle-hole symmetry of the Anderson model combined 
with the Friedel sum rule. The existence of the unitary 
limit can in general serve only as a test of internal numer- 
ical consistency of a particular numerical method and is 
not a direct evidence of a correct description of the Kondo 
physics as argued by THCP. 
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Figure 1: (a) In 2t <C U regime the conductance of THCP (diamonds) compared to results obtained from the two-point 
sine conductance formula-" using variational wave function of Gunnarsson and Schonhammer (full line)£ and CPMC method 
(bullets)?" (b) Conductance for 2t S> U obtained with exact Bethe ansatz (dashed-dotted) 1 — and Hartree-Fock (dashed). In 
numerical calculation t = 1.04(7 was used for the case of variational wave function (full line)^ and CPMC (bullets) * 



In summary, the method used by THCP seems to be promising for the calculation of conductance of systems 
where the interaction is very weak. However, in the limit of strong electron-electron interaction it fails to reproduce 
many-body effects as is, e.g., the Kondo physics in transport through a quantum dot. 
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